Direct observation of mode-coupling instability in two-dimensional plasma crystals 
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Dedicated experiments on melting of 2D plasma crystals were carried out. The melting was 
always accompanied by spontaneous growth of the particle kinetic energy, suggesting a universal 
plasma-driven mechanism underlying the process. By measuring three principal dust-lattice (DL) 
wave modes simultaneously, it is unambiguously demonstrated that the melting occurs due to the 
resonance coupling between two of the DL modes. The variation of the wave modes with the exper- 
imental conditions, including the emergence of the resonant (hybrid) branch, reveals exceptionally 
good agreement with the theory of mode-coupling instability. 
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Strongly coupled complex (dusty) plasmas play an im- 
portant role in the existing hierarchy of soft matter states 
Along with complex fluids and granular media as 
prominent examples of "regular" soft matter, complex 
plasmas represent natural model systems which enable 
remarkably broad interdisciplinary research. The char- 
acteristic length- and timescales in such systems are dra- 
matically stretched (e.g., in complex plasmas - to hun- 
dreds of microns and tens of milliseconds, respectively). 
Therefore, numerous generic processes occurring in clas- 
sical fluids or solids can be studied in greatest detail, 
at the most fundamental "atomistic" level [l|-[3|. Espe- 
cially important here are 2D systems, where the complete 
information about all particles can be obtained at each 
moment of time. The current experimental capabilities 
make 2D complex plasmas ideal for comprehensive ex- 
perimental studies of the atomistic dynamics l2|, y, |5| . 

The investigation of atomistic processes that trigger 
the melting and crystallization is of particular inter- 
est. It is well known that the mechanisms resulting in 
the destruction of the long-range crystalline order in 2D 
systems can be very different from those operating in 
3D crystals: The classical Kosterlitz-Thouless-Halperin- 
Nelson- Young (KTHNY) theory of 2D melting [6] pre- 
dicts two consecutive phase transitions (with an inter- 
mediate, so called "hexatic phase" in between). The al- 
ternative theory relates 2D melting with the formation 
of dislocation chains ( "grain boundaries" ) percolating the 
system. These are generic melting mechanisms which can 
operate in very different strongly coupled systems (sl-flij . 

When studying generic classical phenomena occurring 
in regular liquids and solids, the relevance of a model 
system (be it colloidal suspensions, granular media, or 
complex plasmas) becomes crucial. Of course, careful 
analysis is required in the context of a given phenomenon 
(or, class of phenomena), but one can certainly identify 
essential common principles. In particular, the appli- 
cability of the Hamiltonial approach for the analysis of 
atomistic dynamics (for instance, to investigate kinetics 
of phase transitions) is one of these basic principles. This 
implies that non-Hamiltonian (non-conservative) behav- 
ior peculiar to model systems should play minor role at 



the relevant timescales jl3|- 

Non-Hamiltonial behavior of particle ensembles is of- 
ten a manifestation of asymmetry, or nonreciprocity^ in 
the interparticle interactions. Such asymmetry can be 
considered as a result of the system openness: In the 
presence of strong flow of the ambient plasma the screen- 
ing cloud around each charged particle becomes highly 
asymmetric. These clouds, which are usually referred 
to as "plasma wakes" [13, [HI, play the role of an (ex- 
ternal) "third body" in the interparticle interaction and 
hence make it nonreciprocal [TgI, • One should point 
out that such nonreciprocity is quite universal and can 
be very important for the hydrodynamic interactions of 
colloids [l8[ or droplets in conventional fluid flows [l9| . 

The nonreciprocity of pair interaction in complex plas- 
mas provides a very efficient mechanism of converting the 
energy of the flowing ions into the kinetic energy of mi- 
croparticles [H, H, [20|. In particular, this can be associ- 
ated with dust-lattice (DL) wave modes 

The 

kinetic energy can grow due to resonant particle- wake in- 
teractions when different oscillation modes of interacting 
particles are properly synchronized. Ivlev and Morflll [23| 
suggested that the resonance between DL wave modes in 
a 2D plasma crystal can trigger the mode- coupling in- 
stability which, in turn, can cause the crystal melting. 
By comparing the melting experiments with the molec- 
ular dynamic simulations, this hypothesis has received 
strong conflrmation [24]. However, the direct proof that 
the crystalline monolayers can indeed melt due to the 
mode-coupling instability (and even that the instability 
itself can be observed in experiments) was still outstand- 
ing. One should point out that the other known plasma- 
related scenarios of 2D crystal melting either require the 
presence of localized defects - single particles (or even 
sublayers) above or below the monolayer [25|, or they 
are associated with the charge fluctuations and can only 
operate at very low pressures (well below 1 Pa) [26|, [l^ • 
In this Letter we report on dedicated experimental in- 
vestigation of 2D plasma crystal melting. We explored 
a fairly broad range of experimental parameters cho- 
sen well inside the crystalline regime predicted by the 
generic (KTHNY or grain-boundary) melting scenarios 
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FIG. 1: Sketch of the experimental setup. Microparticles are 
confined above the RF electrode and are illuminated with a 
horizontal laser sheet having a Gaussian profile in the verti- 
cal direction (with a standard deviation cr = 75 /xm). The 
monolayer is levitated well below the peak of the laser in- 
tensity, which results in strong intensity variations of the 
scattered light upon the vertical displacement of individual 
particles. The upward/downward displacement corresponds 
to positive/negative intensity variation. Microparticles are 
recorded from top at a speed of 250 frames per second. 



[12] . The melting was always accompanied by a sponta- 
neous growth of the microparticle kinetic energy, which 
clearly suggests a universal plasma-driven mechanism 
underlying the process. By implementing simultaneous 
measurements of the dispersion relations for three prin- 
cipal DL wave modes, we unambiguously demonstrated 
that the melting occurred according to the mechanism 
predicted by Ivlev and Morfill. The variation of the DL 
modes with the experimental conditions, including the 
emergence of the resonant (hybrid) mode, revealed excep- 
tionally good agreement with the theoretical predictions 
of the mode-coupling instability. 

The melting was systematically investigated in exper- 
iments performed with a (modified) GEC chamber, in a 
capacitively coupled rf glow discharge at 13.56 MHz (see 
Fig. d]). The argon pressure was between 0.4 Pa and 1 
Pa and the forward rf power was between 5 W and 20 W 
(which corresponded to rf peak-to-peak voltage between 
175 V and 310 V). A horizontal monolayer was formed by 
levitating melamine-formaldehyde particles in the plasma 
sheath above the lower rf electrode. We used particles of 
two different sizes, 8.77 ± 0.14 /im and 9.15 ± 0.14 jim 
diameter. A high quality of the monolayer (with no de- 
tected particles above or below) was ensured by using a 
side-view camera. The diameter of the obtained crys- 
talline structure was up to 60 mm, depending on the 
number of injected particles (we explored a range from 
1000 to 15,000). 

One of two different discharge parameters was used as 
the control parameter for the melting: We either grad- 
ually decreased the rf power at a constant gas pressure 
(which lowers the resonance frequency of the vertical con- 
finement /v, see Ref. (ill), or we decreased the pressure 
at a constant rf power (which naturally lowers neutral gas 
friction, but also makes the radial confinement stronger 
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FIG. 2: Example of the instability onset in a 2D plasma crys- 
tal. The shown example is for the experiment performed at 
a rf power of 20 W with particles of 8.77 /xm diameter. The 
pressure of argon was ramped down in steps of 0.02 Pa (every 
:^ 300 s), the instability was triggered at 0.90 Pa. The particle 
trajectories (superposition of consecutive video frames during 
0.04 s) are plotted about 50 s after the critical pressure de- 
crease, the upper and middle panels demonstrate the side and 
top views, respectively. The lower panel shows the profile of 
the particle density in the monolayer measured before the in- 
stability onset (within the stripe marked in the middle panel 
by two horizontal dashed lines). In the upper left corner the 
elementary cell of the hexagonal lattice is depicted with the 
frame of reference chosen in this paper, due to the lattice 
symmetry we consider the wave vector kat0°<^<30°. 



and hence the interparticle distance a smaller) [29]. In 
either case, the melting developed in the same way, as il- 
lustrated in Fig. [2] (see also supplemental movie SI j3Q|): 
Upon decrease of the control parameter below a certain 
threshold the instability was triggered in the central part 
of the crystal (i.e., where the particle density was the 
highest) and then expanded outwards. The instability 
always set in simultaneously in the horizontal and ver- 
tical directions. At the initial stage it was accompanied 
by an exponential growth of the particle kinetic energy, 
eventually resulting in the total disruption of the crystal 
in the unstable zone. Upon increase of the control param- 
eter back above the threshold (usually, with a hysteresis) 
the system always recrystallized. It is noteworthy that 
if the rf power was sufficiently high and/or the particle 



density in the monolayer was low enough the crystal re- 
mained stable at practically any pressure (e.g., at a rf 
power of 20 W, particle density of 160 cm~^, and argon 
pressure down to p 0.4 Pa - the lowest value technically 
achievable with our setup). 

The observed behavior fully agrees with the qualitative 
picture predicted by the theory of mode-couphng insta- 



bility [23|, : It can only set in when the out-of-plane 
and in-plane wave modes intersect, which is only possible 
when /v is sufficiently low or /and the particle density in 
the monolayer (oc a~^) is high enough. The intersection 
should be accompanied by the simultaneous excitation 
of the vertical and horizontal oscillations. On the other 
hand, the instability should be suppressed by sufficiently 
high gas friction, i.e., by increasing the gas pressure. 

The chief result of our paper is the ultimate proof that 
the observed melting is indeed triggered by the mode- 
coupling instability. The unique fingerprint characteriz- 
ing the instability 23|, 2^ is the crossing of the longitu- 
dinal in-plane and transverse out-of-plane DL branches, 
which results in the (wake-induced) resonance coupling 
between them and the emergence of an unstable hybrid 
mode operating in a certain proximity of the crossing. 
This was directly proved by measuring the dispersion re- 
lations of all three principal DL modes. For this pur- 
pose, we used the technique by Nunomura et al. [si^ 
of reconstructing the in-plane dispersion relations from 
thermal fluctuations of particle velocities in the horizon- 
tal plane. In addition, we implemented a novel method of 
particle imaging [22| which allowed us to measure simul- 
taneously the vertical fluctuations and hence to retrieve 
also the out-of-plane mode. The essential advantage of 
this technique is that it does not involve any intermediate 
(model) assumptions about the interparticle interactions, 
and therefore yields "true" dispersion relations (or "fluc- 
tuation spectra") which can be obtained for arbitrary 
experimental conditions (that just allow the individual 
particle tracking). 

In Fig. [3] we plotted the fluctuation spectra for the 
out-of-plane mode (upper panel) separately from those 
for the in-plane modes (second- from-the-top panel), be- 
cause the imaging method implemented for the simulta- 
neous tracking of the vertical and horizontal fluctuations 
yields different background intensities [22]. As regards 
the theoretical curves, in each panel we plotted all prin- 
ciple modes that fit the frequency range shown. 

Figure [3] summarizes results of our analysis. In the 
stable regime (left column) the wave modes are well sep- 
arated and the fluctuation intensity is evenly distributed 
over the branches. In contrast, the unstable regime (right 
column) is characterized by the apparent intersection be- 
tween the out-of-plane and (one of) in-plane branches. 
The intersection occurs in the vicinity of the first Bril- 
louin zone boundary, resulting in the anomalous energy 
release (dark-red "hot spots"). 

The qualitative effects and features of the mode- 




FIG. 3: Evolution of the DL dispersion relations near the in- 
stability onset. The dispersions relations ("fluctuation spec- 
tra") were measured in the experiment performed at argon 
pressure of 0.76 Pa with particles of 9.15 /xm diameter, at a 
rf power of 7 W (left column, stable regime) and 5 W (right 
column, the instability onset). The logarithmic color coding 
is used for the "fluctuation intensity" Ae/A/cA/ (energy den- 
sity, for unit mass, per unit wave number and frequency): Yel- 
low ridges represent stable wave modes, the unstable hybrid 
mode are the dark-red "hot spots" ; the upper panel shows the 
measured out-of-plane mode, the second-from-the-top panel 
depicts the pair of in-plane modes. The solid lines in both 
panels are theoretical results [32] (see text for details). The 
wave number k is normalized by the inverse interparticle dis- 
tance (determined from the first peak of the pair correla- 
tion function; the uncertainty of a was ^ 4% due to density 
variations, see e.g., the lower panel of Fig. [2])- The shown 
range of k corresponds to the first two Brillouin zones at 
— 30° (branches are symmetric with respect to the bor- 
der between the two zones). Below the fluctuation spectra 
the magnified hybrid mode in the unstable regime (marked 
by a rectangular box) is shown. The lower panel illustrates 
the distribution of the "fluctuation energy" Ae/Ak at fixed 
frequency: In the stable case the frequency 12.5 Hz) is 
chosen in the middle between the longitudinal in-plane branch 
and the out-of-plane branch (measured at the border between 
the Brillouin zones), for the instability onset the frequency 
(^11 Hz) corresponds to the hybrid mode. 



coupling instability mentioned above were originally pre- 
dicted by the theory for a ID string isl,!!^, and they all 
remain unaltered also in 2D case. For rigorous quantita- 
tive comparison with the experiments we employed the 
theoretical results for 2D crystals obtained recently by 
Zhdanov et al. 321. The theoretical branches are shown 
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in Fig. [3] by solid lines: The upper curve represents the 
optical out-of-plane mode, the middle and lower curves 
are for the pair of acoustic in-plane modes (longitudinal 
and transverse, respectively). These curves were calcu- 
lated for the normalized dipole moment of the wake q5 
(see Ref. [32|) equal to 0.1. We assumed the Yukawa 
interparticle interaction, with the parameters (viz., the 
particle charge and the effective screening length) de- 
duced from the low-k part of in-plane fluctuation spectra 
(using the method of Ref. [31]). 

Figure [3] shows that the theory yields excellent agree- 
ment with the experiment, both for the stable and un- 
stable regimes. Especially important is that in the un- 
stable regime (right column) the predicted unstable hy- 
brid mode perfectly coincides with the locations of "hot 
spots" in the measured fluctuation spectra. Even subtle 
effects, such as the splitting of "hot spots" at the border 
between the Brillouin zones is very well resolved (see the 
distribution of the fluctuation energy in the lower panel, 
and also the magnified hybrid mode) [33^. 

The fluctuation spectra shown in Fig. [3] were plotted 
for the wave vector k at 6> = 30° (see Fig. [2]), because the 
theory predicts this orientation to be the most unstable 
(see Fig. 4 of Ref. [32]). In order to further highlight 
excellent comprehensive agreement between experiments 
and theory, we plotted additional spectra for different 
angles (see supplemental figure S2 [30|). They clearly 
show that at 6> = 30° the fluctuation intensity of the 
"hot spots" is indeed the highest. 

Note that rapid energy growth after the instability on- 
set results in various nonlinear phenomena which can be 
easily identified in the right column of Fig. [3l The most 
prominent effects are the energy cascades from the "hot 
spots" along the /c-axis (bright red-yellow stripe in the 
second-from-the-top panel, see also the energy profile in 
the lower panel), and the generation of secondary har- 
monics at the double frequency of the hybrid mode (yel- 
low bands at 22 Hz in the fluctuation spectra). 

In conclusion, one should point out that any instability 
operating in 2D plasma crystals and associated with col- 
lective processes [34] should be accompanied by anoma- 
lies in the fluctuation spectra. Therefore, the method 
of simultaneous measurement of all three principal DL 
modes promises to be very fruitful for further investiga- 
tion of the melting scenarios. Especially important would 
be a systematic search for other possible resonant pro- 
cesses that can extend or alter the mechanism reported 
in this Letter. 
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